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1. Introduction 



The dynamics of low energy hadrons are governed by the non-trivial 
structure of the QCD vacuum and the resulting breaking of chiral symmetry. 
Long-range correlations between quarks and antiquarks form condensates, 
whose scale determines constituent quark masses and so the masses of all 
light hadrons. Though the scale has been recently determined for the uu 
and dd through experiments involving tttt interactions, confirming the ~ 
— (235 MeV) 3 anticipated from phenomenology [Ij, we are yet to determine 
a value of the condensate for the not-so-light quarks with any certainty. 

The interest in the value of such a condensate arises in the context of 
QCD sum-rules where the Operator Product Expansion (OPE) is used to 
approximate the short distance behaviour of QCD. In studying currents 
TjjjV (75)^, with qi = s and qj = u, d, the VEVs of uu, dd and ss operators 
naturally arise [21 El II]- The qq condensate for the u and d are expected to 
be close to that in the chiral limit, with even an error of 10% unimportant 
in previous calculations. However, if one considers strange quarks we are 
left with the estimate of Shifman et al [6] that the ss condensate be 
(0.8 ± 0.3) of the uu and dd values. It is the greater precision brought 
about by the studies of Refs. [21 E], for instance, that motivate the need 
to learn about how the qq condensate depends on the current quark mass. 
Our goal here is to illustrate a method for determining this dependence. 



Our aim is to calculate the mass function of the quark propagator for a 
range of current masses. The starting point is the renormalized Schwinger- 
Dyson equation for the quark propagator as depicted in Fig. [U 



2. Schwinger Dyson Equations 
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Fig. 1. Schwinger-Dyson equation for the quark propagator 
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In the Landau gauge we note that Z\ = 1. The inverse propagator S F (p) 
is specified by two scalar functions A and M: 

SfHp) = A(p 2 ,^)^+ M(p 2 )) , (2) 

where the quark mass function M(p 2 ) is renormalisation group invariant. 



2.1. Maris-Tandy Model 

To solve Eq. ([I]) we employ some suitable ansatz for the coupling and 
interaction in Eq. (pQ) which has sufficient integrated strength in the infrared 
to achieve dynamical mass generation. Following Maris et al. [9j [10], we 
employ an ansatz for g 2 D^ v (p — k) shown to be consistent with studies 
of bound state mesons. The renormalisation scheme is one of modified 
momentum subtraction at point taken to be 19 GeV to compare with 
earlier studies [11]. We later evolve \i to the more common 2 GeV scale in 
the MS-scheme. We use: 

- 2 Z 



9 ^ 2 n ( \ 



a 



{l 2 ) D$(q) 



(3) 



where the coupling is described by [HI [TO]: 
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with m t = 0.5 GeV, r = e 2 - 1, 7 m = 12/(33 - 2JV/) and A QC d = 0.234 
GeV. We choose u and to be consistent with meson observables; a typical 
set is u = 0.4 GeV, D = 0.933 GeV 2 . Solutions are obtained by solving for 
A and M of Eq. ([2]), which we may write symbolically as: 

A(p 2 ,v) = Z 2 ( M ,A)-E D (p,A) , 

(5) 

M(;pV(p 2 ,/x) = Z 2 (^,A)Z m m i? ( / u) + S 5 (p,A) . 

The Eg and Sx> correspond to the scalar and spinor projections of the 
integral in Eq. (pQ). For massive quarks we obtain the solution M (later 
called M + ) by eliminating the renormalisation factors Z m via: 

1 + S D (M,A) , 

1 E S (m,A) 



^ 2 (/i,A) 
Z m ((jt,A) 



(6) 



Z 2 (n,A) Z 2 (n,A)m R (fi) ' 

The momentum dependence for different values of mjj are shown in Fig. [2j 
Our purpose is to define the value of the qq condensate for each of these. 
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3. Extracting the Condensate 

At very large momenta the tail of the mass function is described by the 



operator product expansion so that 



M(p 2 
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log (p 2 /A x 2 
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-log(p 2 /A 



7m — 1 



(7) 



where the first term is related to the explicit mass in the Lagrangian, 
by some renormalisation factors. The second term gives the lowest dimen- 
sion vacuum condensate, where C = — (qq). If we included the expression 
to all orders then the scales Ai and A2 would both be equal to Aqcd- 
However, higher contributions to the leading order forms in Eq. ([7]) are dif- 
ferently suppressed, and so Ai and A2 are in practice different. For large 
masses only the first piece is relevant, whereas for m = only the second 
term appears. In this latter case, owing to the power suppression of higher 
orders, A2 is readily determined and found to be equal to Aqcd- Thus in 
the chiral limit we can easily extract the renormalisation point independent 
condensate, C = —(qq), from the asymptotics. 

For non-zero current masses, one can attempt to fit both terms of the 
OPE in Eq. (J7|) to the tail of the mass function, M of Fig. [2j Comparing 
the full mass with m q ^ with that in the chiral limit, one sees how very 
small the contribution of the condensate to the tail is. So while a value for 
the condensate can be extracted, this procedure is not at all reliable. 




Fig. 2. Euclidean mass functions for different current masses, specified at fx = 19 
GeV as labelled. The plot illustrates how on a log-log plot the behaviour dramat- 
ically changes between a current mass of and 3 MeV. 
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Strictly in the chiral limit, we may also extract the condensate using: 
- (qq)^ = Z 2 (/i, A) Z m Gu, A) N c tr D S F (k, fi) , (8) 

with (qq)^ the renormalisation point dependent quark condensate. At one- 
loop, this is related to the renormalisation point independent one by: 

G> iog A^~y m ■ (9) 



which we compare with the asymptotic extraction to good agreement. 

However, for any non-zero quark masses, we cannot apply Eq. (|8|), since 
it acquires a quadratic divergence, cf. Eq. (|7|). Indeed, it is the elimination 
of this divergence that inspired the definition proposed by Chang et al. [TT] . 
which is unfortunately not equal to the condensate of the physical mass 
function and is therefore ambiguous. Consequently, we need a different 
definition, one close to the OPE, Eq. (JT]). 



3.1. Multiple Solutions 

The SDEs, Eq. (pQ), can have multiple solutions being as they are non- 
linear integral equations. In the chiral limit, there exist three solutions 
for Sf(p) and its mass function M(p 2 ). These correspond to the Wigner 
mode , and two non-perturbative solutions of equal magnitude generated by 
the dynamical breaking of chiral symmetry. These we denote by M^' W (p), 
where: 

r M w (p 2 ) = 
M{p 2 ) = { . (10) 

[ M ± (p 2 ) = ±M°(p 2 ) 

Analogous solutions to these exist as we move away from the chiral limit, 
with M~' W restricted to some critical domain < m < m cr . 

Instead of one single solution, we now have three solutions to the same 
model, each with identical running of the current-quark mass (the first term 
in Eq. 0) in the UV region and differing only by their values of the conden- 
sate. Thus, for rnn(fi) < m cr , it is possible to fit Eq. © simultaneously to 
the three mass functions M^, M w . The scales Ai and A2 are found to be 
the same for any current mass within the given model. A2 is equal to Aqqd, 
while Ai is roughly twice as big and depends upon the model parameters. 
The condensates and C w are then determined in an accurate and stable 
way. The results are shown in Figs. |3(a)[ |3(b)| for N f = 0, 4. The error 



bars reflect the accuracy with which the mass functions represented by two 
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Fig. 3. Condensate extracted through simultaneous fitting of the three solutions to 
the fermion mass- function in the Maris- Tandy model with lu = 0.4 GeV and (a) 
Nf = 0, (b) Nf — 4. The functions of the current quark mass are evolved to 2 
GeV in a momentum subtraction scheme. 



terms in the OPE expression, Eq. ([7]), are separable with the anomalous 
dimensions specified. 

We see that within errors the condensate is found to increase with quark 
mass. This rise at small masses was anticipated by Novikov et al. [12] 
combining a perturbative chiral expansion with QCD sum-rule arguments. 
That the chiral logs relevant at very small m q are barely seen is due to the 
quenching of the gluon and the rainbow approximation of Eq. (pQ) . 

We see in Figs. |3(a)| |3(b)| that the M~ and M w solutions bifurcate 
below m cr ~ 43.4(44.0) MeV with uj = 0.4 GeV for Nf = 0(4) respectively. 
But what about the value of the condensate for the physical solution M + 
beyond the region where and M w exist, i.e. mji(fi) > m cr l Having 
accurately determined the scales Ai and A2 in the OPE of Eq. ([7]) in the 
region where all 3 solutions exist, we could just continue to use the same 
values in fitting the physical M + solution alone and find its condensate. 
However, this would make it difficult to produce realistic errors as the quark 
mass increases. 

We see in Figs. 3(a) , |3(b)| that the M~ and M w solutions bifurcate 



below m cr ~ 43.4(44.0) MeV with uj = 0.4 GeV for Nf = 0(4) respectively. 
But what about the value of the condensate for the physical solution M + 
beyond the region where M~ and M w exist, i.e. mji(fi) > m cr ? Having 
accurately determined the scales Ai and A2 in the OPE of Eq. (|7|) in the 
region where all 3 solutions exist, we could just continue to use the same 
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values in fitting the physical M + solution alone and find its condensate. 
However, this would make it difficult to produce realistic errors as the quark 
mass increases. 

Having relaxed the condition that solutions be positive definite, we in 
fact find there exist noded solutions, which have also been discovered recently 
in the context of a simple Yukawa theory [13] . We illustrate this within the 



Maris- Tandy model, for instance with Nf = 4 and uj = 0.4, in Fig. 4(a) We 



note that this noded solution is not limited to the same domain that restricts 
M~ and M w . Thus there exists a solution with a well-defined UV running 
of the quark mass exactly as the M + solution, as far as m = 66.3 MeV. 
While at small quark masses we have all four solutions, at larger masses 
there are still two. Consequently, we can accurately fix the scales Ai and 
A2 of Eq. © at each mn{n), and so determine the condensates as shown 
in Fig. 4(b)[ Indeed, these fits confirm that Ai and A2 are independent of 
mn(fj,). We can then fit the remaining M + solutions shown in Fig. [2] to 
give the physical condensate shown in Fig. |4(b)| for acceptable values of uj 
as determined by [14J. In Fig. [5] we have in fact scaled the quark mass from 
[i = 2 GeV in the (quark-gluon) MOM scheme by one loop running to the 
MS scheme at 2 GeV using the relationship between Amom and A-jg^ for 
4 flavours deduced by Celmaster and Gonsalves [15]. In this latter scheme 
the strange quark mass is ~ 95 MeV |16j . 
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Fig. 4. (a) Momentum dependence of the 4 solutions for the fermion mass-function 
in the Maris-Tandy model with m = 20 MeV at fi = 19 GeV, and (b) Current 
quark mass dependence of the condensates, including the noded solution. We use 
Nf = 4, uj = 0.4 GeV as our parameters. 
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Fig. 5. Condensate for the Maris- Tandy Model with N f = 4, w = 0.4, 0.45, 0.5 GeV 
as a function of current quark mass defined at 2 GeV in MS scheme. 



4. Conclusions 

For the range of parameters considered in the Maris- Tandy modelling 
of strong coupling QCD, we find that the ratio of the condensates at the 
strange quark mass to the chiral limit to be: 

(^) m( MS) =95 Mev/(^) m =0 = (l- 1 ± 0-2) 3 , (11) 

in a world with 4 independent flavours. 

What we have shown here is that there is a robust method of determining 
the value of the qq condensate beyond the chiral limit based on the Operator 
Product Expansion. Of course, as the quark mass increases the contribution 
of the condensate to the behaviour of the mass function, Fig. [21 becomes 
relatively less important and so the errors on the extraction of the physical 
condensate increases considerably. Nevertheless, the method is reliable up 
to and beyond the strange quark mass. Alternative definitions are not. 
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